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Abstract 

We discuss the angular distribution of the lepton pair in the Drell-Yan process, 
hadron + hadron — > 7* + X —* l + + l~ + X. This process gives information 
on the spin-density matrix p( q & of the annihilating quark-antiquark pair in 
q + q — > / + + l~. There is strong experimental evidence that even for unpo- 
larised initial hadrons nontrivial, and therefore the quark-antiquark 

system is polarised. We discuss the possibilities of a general -which 
could be entangled- and a factorising pW>9). "We argue that instantons may 
lead to a nontrivial p( q >V of the type indicated by experiments. 
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1 Introduction 



In this note we discuss the question of factorisation in the Drell-Yan process [1] 

^l + (q + ) + l-(q.). (1) 

Here, h\ and h 2 are the initial hadrons, 7* is the virtual photon, / + , l~ are the final 
state leptons (I = e, fi) and X stands for the hadronic final state particles. The four- 
momenta are indicated in brackets. The basic underlying process is the annihilation 
of a quark-antiquark pair 

q(h) + q(k 2 ) - j*(k) -> l + (q + ) + Z" (?_) , (2) 

which is sketched in Fig. [0 Here we focus on the discussion of reaction (J2J) which is 
the lowest order process in the framework of the QCD improved parton model, see 
for instance [2]. For massless quarks we find that in (J2J) a lefthanded quark q L can 
only annihilate with a righthanded antiquark q R and vice versa. 

The diagram of Fig. ^ is calculated by first evaluating the amplitude for (J2J) and 
folding it then with the parton distributions of the hadrons h±,h 2 . In early theo- 
retical work, one usually assumed that for unpolarised hadrons hi,h 2 the "parton 
beams" delivered by them are also unpolarised. We will call this the no-polarisation 
assumption. In the simplest approximation one furthermore assumes the partons 
q and q to be strictly collinear with the hadrons h\ 2 . This leads to a well known 
angular distribution of the lepton pair in the rest frame of the virtual photon 7*. 
With the polar and azimuthal angles of the outgoing l + , 9 and </>, one gets 

-^ = 7^(1 + 00^). (3) 

O dil 107T ' 

Here and in Eq. (J3J below we use the Collins-Soper reference frame [3] where the 7* 
is at rest and the basis vectors e\ j3 are defined by e 13 = (p 1 ± p 2 )/|Pi ^ P2L w hh 
Pi = Pi/\Pi\- The Collins-Soper frame is obtained from the hi, h 2 cm. system by a 
rotation free boost. 

In general, the angular distribution of the l + is described by three functions 
A, /1, v which may depend on the kinematic variables of (0) 

Ida = fi + ^ cos 2q + n s in 29 cos0+ - sin 2 6> cos20 N ) . (4) 

a as l 47r A + 3 v 2 / 

The LO pQCD result (J3J) implies for the functions A = 1, \i — v — 0. Higher 
order corrections in a s change these values. But within the standard framework and 
using the no-polarisation ansatz one still finds [4] at NLO one relation among the 
coefficients in (JU): 

l-A-2i/ = 0. (5) 
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Figure 1: The generic Drell-Yan process. 



This Lam- Tung relation is almost unchanged at NNLO [5], and even holds for the 
inclusion of parton transverse momentum and soft gluon effects [6]. However, the 
relation (jSJ) is drastically violated in experiments [7-9]. 

In Refs. [10] and [11] two at first sight quite different ideas have been proposed 
giving possible explanations for the violation of the Lam- Tung relation (j5J). It is 
the purpose of the present article to give a short review and a detailed comparison 
of these two approaches. We also sketch a calculation of instanton effects for the 
Drell-Yan reaction and pose some general questions concerning parton factorisation 
and entanglement. 

2 Spin effects and factorisation 

In [12] it was argued on general grounds that the assumption of unpolarised par- 
ton beams from a reaction with unpolarised initial hadrons is questionable due to 
possible vacuum effects. In particular, it was speculated that the fluctuating chro- 
momagnetic vacuum fields which are due to the nonperturbative vacuum structure 
in QCD might lead to a correlated spin orientation of q and q in (j2J) before the an- 
nihilation. This would be in analogy to the Sokolov-Ternov effect [13], well known 
from e + e~ storage rings. 

In [10] this idea was expanded upon and confronted with experiments. A general 
two-particle spin-density matrix for the qq pair in (j2J) prior to the annihilation was 
assumed 



The quantities F i} Gi and Hij are real functions of the invariants of the problem. 



(> 




-{l(8)l+Fy(<r-e*)® 1 
+Gj 1 <g> (tr ■ e*) + {a ■ e*) <g> (tr ■ e*)} . 



(6) 
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Here we work in the qq cm. system and set, 

*i 



63 ~~ \KV 

e * (P1+P2) x g| 

1 " |(Pi* + P 2 *)xe 3 *|' 

e* = e*xe*. (7) 

Such a spin matrix will certainly affect the 7*-production cross section from a qq 
state. The related production matrix in the qq cm. system (the angular distribution 
of / + arises from the contraction with the lepton-production matrix) is, 

r£(fcj, w) = Yl (i:\r\q(K,a,A)q(-kl,f3,B)) 

colours, spins (8) 

x | 5 AA ,5 BB , P%%p,(l*\T\q(kl a', A') q(-k* 17 /?, £')>* ■ 

Here a,/3,a',/3' are the spin indices, A, B, A', B' the colour indices and we have 
assumed no polarisation in colour space. In LO the amplitude for 7* production 
reads, 

W\T\q(k U 0L,A)q(h2,P,B)) = (0|e J^0)\q(h, a, A) q(k 2 ,P,B)) 

= eQ q 5 A BVp(k 2 ) r y fJ ,u a (k 1 ) . 

Here e is the hadronic part of the electromagnetic current. The conventions for 
the Dirac spinors are as in [14] with a = ±1/2 (f3 = ±1/2) representing the quark 
(antiquark) with spin orientation in the direction ±e 3 *. 

With the standard no-polarisation assumption in spin space, one sets 

p (M) |naive = |(l ® 1) = \(S aa ' ■ (10) 

It was shown in [10] that a nonzero correlation coefficient 

H 2 2 — H\\ 



1 + H. 
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could easily explain the experimentally observed deviation from the Lam- Tung re- 
lation. Indeed, defining 

« = -i(l -A- 2v), (12) 
one finds instead of (0) with (jUJ) and (|TT|) 

ft w (k) . (13) 

Here the average is over the parton longitudinal and transverse momenta, see [10]. In 
Eq. (|13|) we do not have an equality sign since higher order perturbative contributions 
give already a (very) small contribution to ft even for k = 0. 
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A good fit to the data of [7, 8] could be obtained with the simple ansatz 

|fc T | 4 

« = |, I a t, K = 0.17, m T = 1.5GeV, (14) 

\kt\ + m T 

where kx is the 7* transverse momentum in the hadronic cm. system. 
In [10] simplifying assumptions for the density matrix (JHJ) were made 

F 2 = F 3 = G 2 = G 3 = H 12 = Hi 3 = H 2 i = H 3 i = . (15) 

This is irrelevant for the Drell-Yan process (J2J). One can easily show that the l + 
angular distribution is only sensitive to the two parameters, ^33 and k. 

Further discussions of possible QCD vacuum effects for the Drell-Yan and other 
reactions were given in [15,16]. 

The problem of the angular distribution in the Drell-Yan process was attacked 
from a different side in [11]. It was pointed out that there can be nontrivial spin 
and transverse momentum correlations even inside an unpolarised hadron. In the 
notation of Ref. [11] the distribution of quarks (with lightcone momentum fraction x\ 
and transverse momentum kix) inside an unpolarised hadron hi (with momentum pi 
and mass Mi) is given by a correlation function $(xi, k^), parametrised as follows: 

*(*!, k 1T ) = fi(xi, k\ T ) £ + hi(x u k\ T ) , (16) 

where the transverse momentum is w.r.t. the plane spanned by the hadron momenta 
Pi and P2 in the hadronic centre of mass frame. In that frame pi and p 2 are pre- 
dominantly in the lightlike n + and n_ directions, respectively. For details we refer 
to Refs. [11,17]. 

The function is the distribution of transversely polarised quarks with nonzero 
transverse momentum inside an unpolarised hadron. The subscript 1 on fi and 
indicates that these functions contribute at leading twist and should not be confused 
with the hadron label. The antiquark correlation function $(x2, k 2 r) is parametrised 
accordingly: 

¥(z 2 , k 2T ) = fi(x 2l k 2 2T ) + hi(x 2 , k\ T ) . (17) 

The quark spin-density matrix can be obtained by projecting $(xi, fcir) onto 
the basis (i[)+r, i^+l), i-e. the right and left chirality components of the good field 
ip + = |7~7 + V ; (see for instance Refs. [18,19]); and analogously for $(£ 2 ,&2t) and 
the antiquark spin-density matrix For given kix and k 2 T one can boost to the 
frame ((7j), which leads (after appropriate normalisation) to 

i{l + F,(<r-e;)} , 
-{1 + G^-e*)} . (18) 
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For simplicity we have suppressed the arguments of the functions. 
From Eq. (fTHj) we arrive at F 3 = = G3 and for i = 1,2: 



fiMS 1 2 ' AM/ 1 11 

G^Jf-^-ei, G^-f^-e;. (19) 

One observes that the function h\ enters in the off-diagonal elements of and 
thus corresponds to RL and LR density-matrix elements. 

In the approach followed in Ref. [11], the qq spin-density matrix is given by the 
tensor product of these two nontrivial one-particle spin-density matrices, 

p M = p (i) g, p ® . (20) 

Clearly, nonzero hj; implies that the standard no-polarisation ansatz ()10|) does not 
hold. Comparison of Eqs. (j2*U|l and (fTHj) with Eq. ©, shows that here Hij = FiGj 
for i,j = 1, 2, 3, and hence H i3 = = H 3i (due to F 3 = = G%). 
The factorisation of the spin-density matrix 

p m i s 

usually implicitly as- 
sumed once factorisation of the dependences on hard and soft energy scales is demon- 
strated for a process. See for instance Ref. [20] for a discussion of factorisation of 
the spin-density matrix in the polarised Drell-Yan process (cf. in particular its Eq. 
(14)). For earlier discussions of issues concerning factorisation for processes where 
transverse momenta play a role see e.g. Refs. [21]. As said, for unpolarised hadrons 
it is standard to choose F t = = G{. Using instead Fi and Gi of Eq. (|19|) in a tree 
level calculation of the Drell-Yan process leads to A = 1,/x = and v ^ 0. The 
general expression for v in terms of h\ is given in Ref. [11], but here we will restrict 
to the case of Gaussian transverse-momentum dependence for illustration purposes. 
We assume that all transverse-momentum-dependent functions are of the form 

/fa, kf T ) = f( Xl ) exp {-R 2 kf T ) — . (21) 

7T 

Moreover, we assume that the width of the Gaussian is the same for f\ and f\ (the 
width will be called R"j) and similarly for h± and (the width will be called R 2 h 
and should be larger than R"j in order to satisfy a positivity bound). This then leads 
to 



K 



v R\ k 2 T * n .-> /.-«/ 



4Rj M X M 2 



exn ( \tf k A ^ l ^M^M (02) 

exp ( - [R h - R f ] —J imXi)fKx2) , (22) 



where R = —(1 — A — 2z/)/4 and e a = eQ q , see Eqs. (|T2jl and (jHJ). We find again 
that the deviation from the Lam- Tung relation arises from an average k (albeit in 
addition to higher order perturbative corrections). In Eq. (|2*2*|l e a denotes the charge 
of the quark with flavour a; the sum is over flavours and antiflavours (indicated by 
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a); and, we have used that f a = f a , i.e. the distribution of antiquarks of flavour 
a inside a hadron h is equal to the distribution of quarks of flavour a inside an 
anti-hadron h. 

Setting Rj = 1 GeV~ 2 and fitting the NA10 data [8] as done in Ref. [11], leads 
to a good fit for R\ - R) = 0.17 Rj and 

J2 a e 2 a hi a (x 1 )hf a (x 2 ) \_ 
E a el /f(si)/f(* 2 ) /"°-° ' [ S) 

where we consider the average over x\ and x%. Assuming n-quark dominance and 
h\ j f\ ~ hi/ f\, this leads to the reasonable result that on average is approxi- 
mately 14% of the size of fi. 

This result is of course dependent on the assumptions, but it serves the purpose 
of illustrating that the data can in principle be explained by a nonzero hj;. Hence, 
in order to experimentally discriminate between the two approaches of Refs. [10, 
11], more data is clearly needed, either from other kinematic regions or from other 
processes. In the next section we will elaborate on what is required and what are 
the opportunities for distinguishing between the two approaches. 



3 Comparison of the two approaches 

In this section we compare the approaches of [10] and [11]. Let us first of all empha- 
sise that the ansatz of [11], given by Eqs. (fTHj) to (j2Dj) is perfectly compatible with 
the general ansatz (jUJ) put forward in [10], but restricts p( q '^ to be factorising. 

There is a further restriction in the ansatz (fTSj) to (|20j) . It requires F 3 = G 3 = 0. 
This comes about as follows. The correlation function <£> in (fTBj) for the hadron hi is 
supposed to depend only on the momenta pi,p 2 , ki (actually only on the direction 
of p 2 ), the correlation function $ for h 2 only on P2,Pi, k 2 . From three four- vectors 
we can only form one axial vector in each case, 

ai 2 = e^p lu p 2p k 1 , 2a . (24) 

To form a pseudoscalar invariant we need all four independent four-vectors 

A = e^p^p 2u k lp k 2a . (25) 

Now _F 3 and G3 are in essence measuring the degree of longitudinal polarisation 
of the quark q and antiquark q. Therefore, due to parity invariance of the strong 
interaction, F 3 and G3 must be pseudoscalar quantities and thus linear in A (|2*3j) . 
But in the ansatz (fTT?|) . (|2T)|) F3 arises from the correlation function $ of hadron hi, 
see fllfij) . and can thus only depend on three four- vectors, from which we cannot form 
a pseudoscalar invariant. Thus, with the ansatz ()18j) to ()20j). F3 must be zero. The 
same holds for G%. Therefore, the ansatz of [11] implies F3 = G3 = and due to 
the factorisation of the qq matrix ^33 = 0. In the general ansatz of [10] the density 
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matrix can from the outset depend on all four four- vectors of the problem, there is a 
pseudoscalar invariant (|25|) available, and _F 3 , G 3 do not have to vanish. Obviously, 
also H33 does not need to be zero in the general approach. 

As mentioned in section 2, the Drell-Yan reaction (j2J) is only sensitive to the 
density-matrix element H33 and the combination k (see Eq. fjlljl ). Therefore, one 
way to check if the restricted form (|2*U|) of is actually realised would be to 
measure if 33. But, as already mentioned in [10], the normalised angular distribution 
(HJ) of the lepton pair is practically only sensitive to k. A factor of 1 + H 33 enters 
in the cross section formula but influences mainly the absolute normalisation. This 
latter effect is difficult to measure due to uncertainties in the quark and antiquark 
distributions and in higher order contributions giving rise to the so-called K-factors. 
Thus we are left with one relevant parameter k. 

Different physical mechanism were proposed in [10] and [11] to produce a non- 
trivial qq density matrix with n ^ 0. In [10] it was suggested that effects of the 
nontrivial QCD vacuum may be responsible for k^O. In [22,23] model calculations 
using the general framework of [11] were performed showing that initial-state gluon 
exchange can produce k^O. 

Let us see if on general grounds we can expect different behaviour for the observ- 
able quantity R (112)1 from these two physical pictures. One possibility for compari- 
son would be to study R as a function of fey- The ansatz given in Eq. (|T1)) -taken 
literally- implies that R ~ k,q for large fey. This is a very different behaviour than 
that expected from an underlying h± function, which is assumed to vanish for large 
quark transverse momentum, in accordance with the ansatz of factorisation of hard 
and soft energy scales in the process. This forces R to vanish (at least at tree level) 
in the limit of large | fey | - Higher order a s corrections may modify this conclusion. 
However, as mentioned the NNLO corrections were shown to be small [5,10]. Their 
(negative) contribution to R was found to be well below 1% for \kr\ values up to 3 
GeV (see figure 6 of [10]). Therefore, one expects R (possibly corrected for the small 
higher order perturbative contributions) to decrease. A constant R, that is both 
positive and large, for large kj- would therefore be irreconcilable with the approach 
of [11]. In the general framework of [10] such a behaviour for R would be possible 
but is certainly not required. 

The dependence of R on the other scale in the process, the lepton pair invariant 
mass (denoted by m 7 * in Ref. [10] and by Q in Ref. [11]), may also be different in 
the two approaches. Unfortunately it is not clear what would be the generic Q 2 
behaviour of R due to vacuum effects. Regarding R arising from nonzero h±, the 
expectation is that it will decrease approximately as 1/Q for large Q. This is based 
on results from Ref. [24], where the influence of soft gluons on similar azimuthal spin 
asymmetries was considered. This means that although R is not power-suppressed 
at tree level, higher order a s contributions effectively give rise to power suppression. 
Note that this is quite different from dynamical higher twist contributions, such 
as discussed in Ref. [25], which typically lead to v ~ 0((k^) /Q 2 ) and therefore, 
are expected to be important only for Q values smaller than the experimentally 
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measured range from 4 GeV up to 12 GeV. In any case, the fall-off or persistence of 
R with increasing Q could be a discriminating feature, similar to the dependence. 

A further possibility to differentiate between the two approaches is to investigate 
a possible flavour dependence of R by varying the types of beams (n ± , p, p). Clearly 
vacuum effects do not favour a flavour dependence. On the other hand, if the ratio 
hj; / fi varies for different flavours and different hadrons, then this could lead to 
an observable flavour dependence. Thus far only tt~ N data have been published, 
although Ref. [9] mentions also to have data for a 7r + beam at the same energy. 

Both vacuum effects and nonzero hi could lead to a R that varies as a func- 
tion of (xi,X2). In addition, observable flavour dependence of this Xi dependence 
would arise if the ratio h^ix^)/ fx(xj) varies differently as a function of x for differ- 
ent flavours and different hadrons. This includes the possibility that h^(xi)/ fi{xj) 
changes its sign as a function of Xj, which would lead to sign changes in R as a 
function of (x\,x<i), even when restricting to only one particular process. 

As a last point in this discussion we mention that since the approach of [11] is 
based on a factorised qq spin-density matrix (j2"Uj) . one can test this type of factori- 
sation by measuring several related processes, such as semi-inclusive deep-inelastic 
lepton-nucleon scattering, where the function enters in combination with other 
functions [26,27]. See also [28,29]. In principle one can determine as many observ- 
ables as unknown functions in order to extract and test the consistency of the 
factorised approach. Needless to say, this is quite a formidable task, but outlines of 
such a scheme have been discussed in some detail in Ref. [17]. 

Clearly it will also be very interesting to compare the predictions of the ap- 
proaches of [10] and [11] for other Drell-Yan type processes, for instance Z produc- 
tion as well as 7* + jet and Z+ jet production in hadron-hadron collisions. 

4 Instanton model 

Thus far we have explained in a rather general way, how a violation of the Lam- 
Tung relation (JHJ) can arise, if the standard ansatz (|1(J|) does not hold. In both 
approaches [10, 11] the strength of the violation follows directly from a comparison 
with the experiment. It would be of great interest to calculate the relevant parameter 
describing the asymmetry (namely k) in a certain model. For the approach of 
Ref. [11] this has been done in Ref. [22,23] using a spectator model. It was shown 
that initial-state gluon exchange could give rise to a nonzero and a corresponding 
R for the pp and 7T~p initiated Drell-Yan process. In this section we want to take a 
different approach, namely to outline a model calculation that is along the lines of 
Ref. [10]. 

We already discussed the possibility that in a nontrivial vacuum the spins of 
the partons might be correlated. An intriguing possibility to describe the vacuum 
structure is given by instantons. Instantons [30] are nonperturbative fluctuations 
of the gluon fields and well known to induce chirality-violating processes, absent in 
conventional perturbation theory [31]. Especially this feature of instanton- induced 
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processes is a strong motivation to study the role of instantons as a source of spin 
correlations. Along similar lines, various remarkable effects induced by instantons 
were investigated. One can find for instance in Ref. [32] an estimate of certain single 
spin asymmetries and in Ref. [33] an estimate of the Pauli form factor of the quark. 
One also expects an impact of instantons on the question how the proton spin is 
built up from the spin and angular momenta of the constituents, see e.g. [34] or 
the review [35]. Recently, an estimate of an azimuthal spin-asymmetry induced by 
instantons in semi- inclusive deep-inelastic scattering was presented in Ref. [36]. 

In a rather qualitative way, the Drell-Yan process was already investigated in an 
instanton background in [37]. There it was argued that even in the limit of high 
energy, instantons may lead to sizable effects, not suppressed by inverse powers of 
the energy. But it should be mentioned that spin effects do not play any specific 
role in [37]. 

Here, we want to emphasise that instantons might indeed violate the naive ansatz 
(jl(Jj) via some additional terms to the spin matrix. The generic instanton-induced 
chirality- violating process which contributes in the Drell-Yan case reads for rif active 
flavours (see Eq. (j2J) for the similar process in usual perturbation theory), 

q L + q L ->7* + (n/-l) [q R + q R ] + n g g (26) 

and is sketched in Fig. |21 The indices indicate the helicity of the quarks and anti- 
quarks in the process. Of course, the process with R and L exchanged everywhere 
-induced by antiinstantons- must also be taken into account. The important point 
in our approach is not the significant complication of the final state in (J26|) which 
contributes to the final state X in (|I|L but the different helicity structure in the 
initial state. 

We mentioned in section 1 that neglecting quark masses in the process (J2J) only 
a quark and an antiquark with different helicities couple to the photon. So one 
can split the process (|26p into two stages: during the first stage, the quark (or the 
antiquark) will change the helicity and afterwards the quark (antiquark) will interact 
in the usual way with the antiquark (quark). The final state will only change the 
size of the whole instanton contribution but not the structure of the related spin 
matrix p( g '^. 

In Fig. |21 we show the two amplitudes contributing to the process (J26|) and for 
each amplitude the split into two stages. In the left part (labelled in Eq. (I27|) 
with (t)) the incoming quark changes the helicity. The right part (u) is of course 
similar but the incoming antiquark changes the helicity. It is sketched that both 
amplitudes factorise into an instanton part described by the coefficients a* 7 -* and 
and chirality-conserving amplitudes, namely (j*\T\q R q L ) or ("y^\T\q L q R ). 

For the simplest instanton-induced process with rif — 1 and n g = (see [38] for 
a detailed calculation of the related process in lepton-hadron scattering) one would 
expect a trivial connection between and For the general case this will 

change because of the more complex kinematics, related to the additional momenta 
of the final state partons. 
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+ 




= \aW(y;\T\q R q L ) + 6« {^\T\q L q R ) \* 
Figure 2: The instanton- induced process q L + q L — > 7* + (nj — 1) [g H + + n g g . 



The important point is that the two processes shown in Fig. El lead from the 
same initial to the same final states. Therefore these amplitudes must be added 
coherently. This gives in the cross section a term 



fj,LL ii LL 



2 uLL^~ 1 vLL 



,C0l 2 



W /) & (/) *(7;|T|g fl gJ(7:|T|g i ^* + a( / )*6( / )( 7 ;|T|g i ^)( 7 :|T|g,gJ* 

+\b {I) \^i;\T\q L w(i:\T\q L q R y ■ m 

Comparing the general amplitude © with the instanton-induced one (j2*7|) . we get 
the following expressions for the density-matrix elements (the factor 1/4 arises from 
the averaging over the initial state helicities), 



o (I) 

Hrl rl 



,( 7 )| 2 



o {I) 

1 HLRLR 



o (I) 

? FRLLR. 



o {I) 

? Hlrrl 



• (28) 



For the calculation of the spin matrix, we have to add the contribution from the 
usual process without instantons in the background. The naive expectation related 
to (HUJ) is p MM | naivo = pR L R L \ nSlivc = 1/4 and p RLLR \ nSlive = Plr rl I naive = °- Adding 
this to the instanton-induced contribution we get finally, 



P 



in: l r; ~r Plrhl 



2 Re (oW>*) 
'2 + |a( 7 )| 2 + |6( 7 )| 2 ' 



(29) 



An estimate of n in the simplest case where nf = 1 and n g = leads to k 7^ and 
we expect the same to be true in general. 



11 



We want to mention that one can expect contributions to k also from instanton- 
induced processes without any additional partons in the final state. In this case an 
instanton-antiinstanton pair is located only on one side of the cut appearing in the 
contributions to the cross section (in contrast to the squared amplitudes in Fig. 121 
where one instanton will appear on each side of the cut). Hence, the quark and 
the antiquark will change the helicity on one side of the cut and we will also get 
off-diagonal contributions. 

We summarise: The flipping of the helicity of one quark or antiquark in the 
initial state which occurs in the instanton-induced contribution to the Drell-Yan 
process should give rise to a nonzero matrix Hij in (jHJ) and finally to k. As already 
mentioned, the Drell-Yan process (J2J) is not sensitive to Fi and Gi, hence we can not 
say anything about instanton-induced contributions to Fi and Gi. A more careful 
analysis of the instanton-induced contributions to the Drell-Yan process including 
the complete final state in (|26|) is beyond the scope of the present paper. This 
and the question whether the instanton induced processes lead to a factorising or 
entangled qq density matrix will be investigated elsewhere. 

5 Summary 

In this paper we have discussed the angular distribution of the lepton pair in the 
Drell-Yan process (JTJ). We considered the lowest order reaction (j2J) and studied 
the influence of the quark- antiquark spin-density matrix on the lepton's angular 
distribution. 

It is well known that a trivial spin-density matrix (jl(Jj) is disfavoured by experi- 
ment (see [7-9]). Experiments are well described by qq spin-density matrices having 
the coefficient k^O. This can be achieved by a qq density matrix which is factoris- 
ing into nontrivial q and q single-particle density matrices, as assumed in the ansatz 
of [11]. The ansatz of [10] is perfectly compatible with this, but would allow also 
for truly entangled qq pairs, that is a two-particle spin-density matrix which cannot 
be written as a tensor product of one-particle matrices. We have made a detailed 
comparison of the approaches [10] and [11] and have shown how they are related. 
We have discussed the underlying physical ideas and have outlined ways to check 
these ideas experimentally. 

We have discussed instanton effects on the quark-antiquark density matrix and 
argued that these could induce spin correlations of the type indicated by experi- 
ments. The question whether instantons lead to a factorising or an entangled qq 
density matrix will be studied elsewhere. 

We think that it is an important question to follow up how to determine from 
experiments the complete qq density matrix. For this other reactions besides the 
Drell-Yan process clearly are needed. We have given some discussion of this issue in 
section 3. It would be fascinating if the qq density matrix turned out to be entangled. 
Thus, in this article we want to pose the question: can there be entanglement at the 
parton level in hadronic reactions? 
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